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Abstract 

We report on the nonlinear coupling between the mechanical modes of a nanotube resonator. 
The coupling is revealed in a pump-probe experiment where a mode driven by a pump force is 
shown to modify the motion of a second mode measured with a probe force. In a second series 
of experiments, we actuate the resonator with only one oscillating force. Mechanical resonances 
feature exotic lineshapes with reproducible dips, peaks, and jumps when the measured mode is 
commensurate with another mode with a frequency ratio of either 2 or 3. Conventional lineshapes 
are recovered by detuning the frequency ratio using the voltage on a nearby gate electrode. The 
exotic lineshapes are attributed to strong coupling between the mechanical modes. The possibility 
to control the strength of the coupling with the gate voltage holds promise for various experiments, 
such as quantum manipulation, mechanical signal processing, and the study of the quantum-to- 
classical transition. 
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The nonlinear nature of mode coupling lies at the origin of a wide variety of phenom- 
ena p]-[T2]. including mechanical synchronization, mechanically-induced transparency, and 
vibration localization. In the case that two modes operate at well separated frequencies, 
the effect of the coupling is usually modest and the oscillators move essentially in an in- 
dependent manner. The coupling between two modes is expected to become strong when 
the ratio between their resonance frequencies is an integer n [13j. Perturbation theory then 
predicts that the motion of one oscillator strongly affects the motion of the other oscillator, 
and vice-versa^ through nonlinear forces of order n (see below). Moreover, the resonance 
lineshapes are expected to be peculiar [13j. Strong coupling, also called internal resonance, 
has not been observed in nanomechanical resonators thus far, because the ratio between the 
resonance frequencies is usually not an integer. 

Resonators based on carbon nanotubes [T4l-[T6] provide a unique platform to test mechan- 
ics at the nanoscale. A nanotube behaves like a semi- flexible polymer in the sense that it can 
bend and stretch to large extents Consequently, nonlinearities in nanotube resonators 
are important and result in unusual behaviours [ 181420] . In this work, we take advantage 
of the mechanical flexibility of nanotubes to achieve strong coupling between mechanical 
modes. Indeed, the static shape of a nanotube can be deformed to a large extent with the 
voltage applied on a gate electrode. This enables us to tune the resonance frequencies p!4l[2T] 
in order to make two modes commensurate. In addition, the oscillation is easily driven to 
large amplitudes [14j so that nonlinear forces, including coupling forces, are sizeable. 

In this Letter, we study the lineshape of mechanical resonances as a function of gate 
voltage. Lineshapes become exotic, featuring reproducible dips, peaks, and jumps, when the 
measured mode is commensurate with another mode with a frequency ratio of either 2 or 3. 
Conventional lineshapes are recovered by detuning the frequency ratio with the gate voltage. 
These results agree with the predictions of strong coupling. The coupling is attributed to 
motion-induced tension, that is, the oscillation of one mode induces a mechanical tension in 
the resonator that affects the dynamics of the other mode, and vice-versa. 

We employ conventional techniques for the fabrication and the measurements of nanotube 
resonators. Figures 1(a) and (b) show that the nanotube is contacted to two electrodes and 
is suspended over a trench with a gate electrode at the bottom. The nanotube is grown 
by chemical vapour deposition in the last step of the fabrication process in order to reduce 
contamination [22j (supplementary section I). We check with a scanning electron microscope 
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that only one nanotube is suspended over the trench. The mechanical motion is driven and 
detected using the two-source and the frequency-modulation (FM) mixing methods. The 
two-source method which enables a direct measurement of the amplitude of the motion, 
is used to record resonance lineshapes, whereas the FM method [23] is better at detecting 
small signals so we employ it to map resonance frequencies as a function of gate voltage 
Vg (supplementary section II). Measurements are performed between 60 and 70 K to avoid 
Coulomb blockade at low temperature p!5l [T6] . 

The nanotube resonator already begins to exhibit Duffing nonlinearities at low driving 
force Fd [24]. Figures 1(c) and (e) show two resonance lineshapes at the lowest Fd for 
which we obtain a good signal-to-noise ratio. The two resonances correspond to a single 
mode at different values of Vg. The quality factors are 230 and 350 in Fig. 1(c) and (e), 
respectively. Upon doubling Fd^ a hysteresis emerges, marking the onset of the nonlinear 
regime [Fig. 1(d), (f)]. An estimation of the motional amplitude yields values between Inm 
and 9nm in Fig. l(c-f) (supplementary section III). Interestingly, the asymmetry of the 
resonance is different between Fig. 1(d) and (f), which indicates different signs of the Duffing 
force. The sign change occurs around Vg = 1.9 V. 

The resonance frequencies can be tuned with Vg by an amount that is different for each 
mode [Fig. 2(a),(b)]. The resonance frequency variation is attributed to the mechanical 
tension that builds up in the nanotube as it bends towards the backgate upon increasing 
Vg [m [25]. The amount of the variation depends on the shape and the direction of the 
mode. Finite element simulations can qualitatively reproduce the measured Vg dependences 
of the different resonance frequencies [Fig. 2(c)] without any free parameters using the static 
shape of the nanotube imaged with a scanning electron microscope (supplementary section 
VII). These simulations show that the static deformation of the nanotube towards the gate 
electrode is as large as 50 nm for 1^ = 4V [Fig. 2(d)]. For the simpler case of a straight 
nanotube, we can describe the Vg dependences of the resonance frequencies in a satisfactory 
way using the Euler-Bernoulli equation (supplementary section VIII); the static deformation 
is 17 nm for 1/^ = 4 V [Fig. 2(e)]. 

Coupling between the modes can be observed in a pump-probe experiment [5]. Specifi- 
cally, we apply a force at frequency f probe to probe one mode using the FM method. The 
current of the probed mode is continuously monitored while sweeping the frequency fpump 
of a second force [Fig. 3(a), (b)]. The sweep in fpump is repeated for various values of Vg 



[Fig. 3(c)]. The current of the probed mode is found to change when fpump matches the 
resonance frequency (or the harmonic) of another mode [by comparing Fig. 3(c) and sup- 
plementary Fig. S3(a)]. This unambiguously demonstrates that the modes of our nanotube 
resonator are coupled. 

When only one mode is actuated, we observe discontinuities in maps of the resonance 
frequency as a function of Vg [Fig. 4(e), (f)]. The discontinuities are accompanied by exotic 
resonance lineshapes [Fig. 4(g-j)]. These features often occur when the resonance frequency 
of the measured mode is equal to that of another mode multiplied by 2, 3, 1/2, or 1/3 
(as indicated by the black arrows in Fig. 2(b); [26j). Upon detuning the frequency ratio 
using Vg^ conventional resonance lineshapes are recovered [Fig. 4(a), (d)]. Another way to 
retrieve regular lineshapes is to reduce the driving force (supplementary section XIII). We 
also observe exotic lineshapes [indicated by gray arrows in Fig. 2(b)] without being able to 
identify the second mode; we speculate that the second mode is not detectable with the 
mixing technique or that it oscillates in a frequency range that has not been probed. 

These experimental findings are consistent with the theory of strong coupling between 
mechanical modes in a resonator p!3l [27] . The observation that strong coupling occurs for 
a frequency ratio of 2 or 3 implies that quadratic and cubic nonlinear forces are important 
and that the equation of motion for mode i is of the form 



= -cj-Zi - 7— - ^2^^ -asz- - 2_^PjkZjZk- 2_^ejkiZjZkZi + g (1) 

with Zi the motional amplitude, t the time, uJi the angular resonance frequency, and g 
the eflFective force normalized by the mass [13j; 7, a2, 0^3, /3jk^ and cjki are various constants. 
We omit the nonlinear damping force for simplicity [18j. Mode i couples to modes fc, and 
/ through the forces ZjZk and ZjZkZi (supplementary section IX). 

Quadratic and cubic nonlinear forces (zf , zf^ ZjZk and ZjZkZi) naturally emerge from the 
tension in the beam that is induced by motion - the beam is stretched and compressed 
periodically in time because it is clamped at both ends. The zf and zf forces are responsible 
for the hystereses and the asymmetric resonance lineshapes in Fig. 1(d) and (f). The upward 
asymmetry in Fig. 1(d) is associated with the cubic zf force, since motion-induced tension 
leads to a positive coefficient a^. When the static deformation of the beam Zg becomes 
sizeable, the quadratic zf force can lead to a reversal of the asymmetry [28l[29]. We estimate 
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from the asymmetries in Fig. 1(d) and (f) that Zs is 2.8 and 13 nm at = 1.5 and 4V, 
respectively (supplementary section XI). This is in fair agreement with the calculation in 
Fig. 2(e), which supports that the nonlinear zf and forces originate from motion-induced 
tension. We estimate that these forces are 3 orders of magnitude larger than electrostatic 
nonlinear forces [28] and thus neglect the latter (supplementary section XII). The coupling 
forces ZjZk and ZjZkZi are intimately related to the z'f and zf forces, since they all arise in the 
same way from the Euler-Bernoulli equation (supplementary section IX). It is thus likely that 
the modal coupling in our experiment is also due to motion-induced tension. In other words, 
the coupling is mediated by the tension generated by the oscillation of one mode, which 
affects the dynamics of the other mode, and vice-versa. The solutions of the equations of 
motion that describe motion-induced tension [Eq. (1)] are characterized by exotic lineshapes 
for the case of commensurable resonance frequencies [13j. The lineshapes are sensitive to 
the coefficients of the coupling forces in a critical fashion. A detailed comparison between 
experiment and theory is not possible at the moment, since the coefficients depend on the 
static shape of the nanotube, which is not known precisely enough. 

The exotic lineshapes in nanotube resonators are analogous to Fermi resonances observed 
in the infrared and Raman spectra of molecules [301 l3T j . When the frequency of a vibrational 
mode of a molecule is twice as large as that of another mode, energy can be transferred from 
one mode to the other. This leads to a mixing of the eigenfunctions and to unusual spectra. 
However, the coupling between the vibrational modes cannot be externally tuned as in 
nanotube resonators. 

The mode coupling force can be made larger in nanotube resonators than in resonators 
made from other materials, since the coupling force scales inversely with the fourth power 
of the resonator length (supplementary section IX) and that nanotube resonators can be as 
short as ^ 100 nm [32l [33j. Mode coupling is further enhanced by the excellent material 
characteristics of nanotubes, since the coupling force is linearly proportional to E/p (sup- 
plementary section IX) and that nanotubes have a high Young modulus E and a low mass 
density p. 

The achievement of strong coupling combined with the possibility to tune its strength 
open up many possibilities. Such coupling may lead to sizeable signatures in the quantum- 
to-classical transition of a mechanical resonator [34j. In the quantum regime, it may allow for 
the manipulation of energy quanta between different mechanical modes using gate voltage 



pulses. Classically, the transfer of energy between mechanical modes could be made faster 
than the energy relaxation time, which is interesting for high-speed signal operation [TTl 
[35] - [37] . The nonlinear nature of strong coupling is expected to give rise to non-intuitive 
behaviours that have not been tested thus far [13j. A striking example is that driving one 
of two coupled modes can cause the second mode to reach a higher amplitude than that of 
the driven one. 

When finalizing the manuscript, we became aware of the paper by Antonio et al. ^38j 
that reports on strong coupling in a 0.5 mm long micromechanical resonator. The frequency 
of the modes is tuned by increasing the driving force (through the Duffing force). The 
possibility to tune the resonance frequencies of a nanotube resonator with a gate voltage is 
more convenient for practical use. 
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FIG. 1: Device characterization, (a) Coloured scanning electron microscopy image of the device 
measured in this work with source (S), drain (D), and gate (G) electrodes. The nanotube position 
is represented by a dashed line and the clamping points are indicated by arrows. The suspended 
length of the nanotube is 1.77 /im and the depth of the trench is 370 nm. Scale bar: 600 nm. 
The image is recorded after the measurements of the resonator, (b) Schematic sideview of the 
device, (c-f) Mechanical resonances for small and large driving forces obtained by measuring the 
mixing current {Imix) cis a function of the driving frequency (/) with the two-source technique. 
The driving force is electrostatic and is proportional to the oscillating voltage V^^ applied to the 
gate electrode. = 0.2mV and Vg = 1.5V in (c); = 0.4mV and Vg = 1.5V in (d); 

= 0.2 mV and = 4V in (e); = 0.4 mV and = 4V in (f). For the detection, we apply 
an oscillating voltage {Vg^ = 0.056 mV) to the source electrode. Black (red) curves correspond to 
upward (downward) sweeps. 




FIG. 2: Tuning resonance frequencies, (a) Map of resonance frequencies as a function of Vg 
(by measuring I^ix as a function of / and Vg for V^^ = 2mV with the FM technique). We 
clearly discern three modes, while a fourth one is weaker and is indicated by arrows. Colour scale: 
(black) to 0.1 nA (red), (b) Schematic of the map of resonance frequencies as a function of 
Vg. The four modes are represented by plain lines and labelled M, O, and P. Dashed lines 
correspond to the resonance frequencies of these modes multiplied by 2, 3, 1/2, or 1/3 (the values 
are indicated in the labels). Black arrows point to regions where lineshapes are exotic and two 
modes are commensurate. Grey arrows point to exotic lineshape regions for which we cannot 
assign the coupled mode, (c) Finite element simulation of the map of resonance frequencies as a 
function of Vg obtained with ANSYS. The dashed hne correspond to a mode that we have not 
detected. The simulations show that this mode has one node (while the others have either zero or 
two nodes) and thus cannot be detected due to symmetry reasons. Inset: schematic of the static 
shape of the nanotube when = V. The deformation in the transverse direction is exaggerated 
with respect to the nanotube length. The largest deformation is ^ 40 nm (supplementary section 
VII). (d) Static displacement of the center of the resonator (zg) calculated with ANSYS using the 
static shape of the nanotube (when = OV) depicted in the inset of (c). (e) Zg calculated from 
the Euler-Bernoulli equation for a straight nanotube (plain line). The dashed line corresponds to 
the result calculated with ANSYS for the same straight nanotube. 
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FIG. 3: Mechanical coupling measured in a pump-probe experiment, (a) Representation of the 
two drive frequencies used in the pump-probe experiment. The pump force is swept in frequency, 
whereas f probe is set to match the resonance frequency of the lowest mode (M). (b) Normalized 
mixing current of the probed mode as a function of fpump at = 3.6 V (measured with the FM 
technique). Before the scan, we set f probe so that the current is maximal (I^i^)- P^^^ 
measured current divided by Z^^^,. The oscillating voltage of the pump is 5.6 mV and the FM 
oscillating voltage of the probe is 2 mV. (c) Normalized mixing current of the probed mode as a 
function of fpump and Vg using the same parameters as in (b). The line graph in (b) is marked 
with a dashed line. Colour scale: Imix — (dark red) to Imix — 1 (white). 
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FIG. 4: Resonances when the measured mode is commensurate or nearly commensurate with 
another mode, (a-d) Resonance hneshapes of one mode for different Vg measured with the two- 
source technique. The driving voltage applied to the gate electrode is V^*^ = 1.1 mV and the voltage 
applied for the detection on the source electrode is V^*^ — 0.56 mV. Black (red) curves correspond to 
upward (downward) sweeps. = 2.4 V in (a); = 2.6 V in (b); = 2.8 V in (c); = 3 V in (d); 
(e) Map of the resonance frequency as a function of Vg (obtained by measuring Imix ^ 

function 

of / and Vg for V^^ = 1.1 mV and V^^ = 0.56 mV using the two-source technique by increasing 
/). Colour scale: (black) to 1 nA (red), (f) Map of the resonance frequency as a function of 
Vg (obtained for V^^ = 1.1 mV and V^^ = 0.28 mV using the two-source technique by decreasing 
/). Colour scale: (black) to 0.7 nA (red), (g-j) Resonance lineshapes for different modes and 
different values of Vg measured with the two-source technique. V^^ = 1.1 mV, V^^ = 0.56 mV, 
and Vg = 1.98V in (g). = 17mV, 1//^ = 1.1 mV, and Vg = 3.41V in (h). = 1.7mV, 

= 0.28 mV, and Vg = 1.86V in (i). = 5.6mV, 1//^ = 0.56mV, and Vg = 3.6V in (g). 
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I. DEVICE FABRICATION 

Our nanoresonators consist of a suspended carbon nanotube clamped between two metal elec- 
trodes, as depicted in Fig. 1(a) and (b) of the main text. The devices are fabricated as follows. A 
trench is etched into a highly resistive Si wafer coated with Si02 and Si3N4. W and Pt are evapo- 
rated into the trench to create a gate electrode (G). In a second lithography step, a continuous line 
is exposed across the trench. After a deposition of W/Pt and hft-off, the hne results in the source 
(S) and drain (D) electrodes separated by the trench (these electrodes are electrically isolated from 
the gate due to the undercut profile of the Si3N4/Si02 substrate). W and Pt are chosen because 
of their high melting points that allow the growth of carbon nanotubes. An island of catalyst is 
patterned on the drain (or source) electrode using electron-beam lithography. Nanotubes are grown 
by chemical vapour deposition from these islands. In about 1 out of 20 cases, a nanotube grows 
across the trench and establishes electrical contact between S and D. This growth is the last step 
of the fabrication process so that nanotubes are not contaminated with residues from resists and 
chemicals [T]-[3]. The device we present in the main text has a length of 1.77 /xm. The separation 
between the nanotube and the gate electrode is 370 nm. The rather large roughness of the S and 
D electrodes in this device does not allow us to measure the nanotube radius with atomic force 
microscopy. 

II. MEASUREMENTS DETAILS 

Our measurements are carried out at pressures typically below 10~^ mbar and temperatures 
between 60 and 70 K. In order to clean the nanotube surface, we perform a current annealing step 
every day (6/xA for 300 s). We observe only very minor variations of the electrical conductance 
and mechanical resonance frequencies of the nanotube from day to day. 

1.^ 



The chip containing the device is mounted on a printed circuit board, dc and ac voltages are 
added through a bias tee outside the chamber. The low frequency mixing current is measured from 
the drain (D) electrode and is low-pass filtered through a capacitor to ground (InF). 

We discuss first the frequency mixing (FM) technique [4J. A driving voltage V^^ is applied to 
the source electrode. Modulating the frequency (with a modulation rate of 671 Hz and a frequency 
deviation of 100 kHz) results in a mixing current {Imix) at 671 Hz. The gate electrode is biased 
with a dc voltage Vg to tune the resonance frequencies. 

In the two-source technique [5j, we apply the driving voltage V^^ to the gate in addition to a 
dc voltage Vg. The motion of the nanotube is detected by applying a second, smaller voltage V^^ 
to the source. The two oscillating voltag es are slightly detuned, and the amplitude signal of Imix 
is measured at the detuning frequency {5uj/2ti — 10 kHz). 

The modulus of I^ix measured with the two-source technique has the form 

1 dG { C \ 

Imix = l^s""^ (^y'^^COs(5cjt - LPe) + ^^yg-^ COs(Sut - (fE - ^m)J (2) 

where G is the conductance of the nanotube, (/^^ is the phase difference between the voltages 
applied to source and gate, t is time, zq is the mechanical amplitude, C is the capacitance between 
the nanotube and the gate, C is its derivation with respect to the nanotube displacement, and 
(fM is the phase difference between the nanotube displacement and the driving force. 

The measurements in the paper of Imix as a function of the drive frequency / give a resonance 
lineshape that is to a rather good approximation proportional to the response of the motional 
amplitude as a function of /, since the purely electrical component of Imix (first term in Eq. [2]) is 
much lower than the mechanical component (second term in Eq. [2]). 

We verify that the harmonics of the RF sources (signals at 2, 3, or 1/2 times the drive frequency) 
can be neglected. These harmonics are far below the smallest driving voltage for which we can 
detect a resonance. Namely, the voltage of the harmonics is typically 1000 times lower than V^^. 



III. ESTIMATION OF DYNAMICAL AMPLITUDE 

Equation [2] allows estimating the motional amplitude of the resonator by comparing the signal 
on resonance, Imax^ to the purely electrical background far from resonance, Iback [5j. Using the 
approximation C = \n(2{d-z)/r) ' ^^^^ 

U 





FIG. 5: Maps of resonance frequencies as a function of gate voltage, (a) Two-source measurement 
at low driving force (obtained by measuring I^ix as a function of / and Vg with V^^ = 1.7 mV and 
V^^ = 0.3 mV). Colour scale: (blue) to InA (red). For comparison, we plot the schematic of the 
modes detected with the FM technique in (b) and the results of the ANSYS simulation in (c). The 
inset in (c) shows the static nanotube shape (measured in the absence of an applied dc voltage) 
that is used in the simulation. 



with d = 370 nm the equilibrium distance between the nanotube and the gate electrode. Since 
we cannot measure the diameter of the nanotube due to the large surface roughness of the electrodes 
in the studied device, we use a typical value for the radius (r = 1.5 nm). We find the following 
values for the maximum mechanical amplitudes zq in Fig. 1 of the main text: zq c:^ 3.2 nm in 
Fig. 1(c), Zq 2^ 8.7 nm in Fig. 1(d), zq 2^ 0.9 nm in Fig. 1(e), and zq 2^ 2.1 nm in Fig. 1(f). 



IV. MAPS OF RESONANCE FREQUENCIES AS A FUNCTION OF GATE 
VOLTAGE: COMPARISON BETWEEN THE TWO-SOURCE AND THE FM 
TECHNIQUES AT LOW DRIVING FORCES 

We carry out measurements of the resonance frequencies as a function of Vg with the two-source 
technique at low driving force [Fig. [5]^a)]. Two modes are detected with a large signal, a third one 
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with a small signal. To facilitate a comparison between the two measurement techniques, we plot 
the schematic of the modes detected with the FM technique in Fig. ^h). We find that the modes 
producing a large signal with the two-source technique are those we label N and P. The mode 
producing a weak signal is identified as O, while M is not detected at all (it does show up at larger 
driving forces). 

The ANSYS simulation helps understanding the relative strengths of the signals [Fig. ^c)]. 
The details of the simulation are discussed in section VII. N and P correspond to modes moving 
essentially in the plane orthogonal to the gate electrode. Since the two-source method measures 
the oscillation of the nanotube-gate capacitance, the signal of modes N and P are expected to be 
large, in agreement with the experiments. M and O are modes moving essentially parallel to the 
gate electrode. As such, the corresponding signals are expected to be small, which also agrees with 
the measurement. 

ANSYS simulations indicate that modes M, N, O, and P have either or 2 nodes. The mode 
with 1 node moving perpendicular to the gate electrode is predicted to appear between N and O 
[dashed line in Fig. [5]^c)]. This mode is not detected because the oscillation of the capacitance is 
(nearly) zero due to the symmetry of the mode shape. ANSYS predicts that the mode with 1 node 
parallel to the gate has a frequency larger than that of modes O and P. 

V. MAPS OF RESONANCE FREQUENCIES AS A FUNCTION OF GATE VOLT- 
AGE: COMPARISON BETWEEN LOW AND HIGH DRIVING FORCES 

In the maps of resonances as a function of Vg the number of detected resonances depends on 
the driving force F^. Figure |6]^a) shows the spectrum obtained with the two-source method for a 
low driving force. Upon increasing by a factor 10, many more resonances appear Fig. [6j(b)]. 
A comparison of the resonance frequencies reveals that the additional resonances are almost all 
harmonics of the four modes M, O, and P (see following section). 

The same scenario develops for measurements with the FM technique at high driving force: 
additional resonances are detected [Fig. [7|^a) and (b)] and identified as harmonics of the four 
modes (see following section). Here, the number of discernible harmonics is even larger than that 
measured with the two-source technique. 




FIG. 6: Map of resonance frequencies as a function of gate voltage, (a) Two-source measurement 
at low driving force (obtained by measuring I^ix as a function of / and Vg with V^^ = 1.7 mV and 
V^^ = 0.3 mV). Colour scale: (blue) to InA (red), (b) Same measurement with a larger driving 
force (V^^ = 17 mV and V^^ = 1.1 mV). White arrows point out the lowest mode which is faintly 
visible. Colour scale: (blue) to 7nA (red), (c) Schematic of the modes detected with the FM 
technique for comparison. 

VI. HARMONICS 

In Fig. [8} we plot the four modes and their harmonics on top of each other by dividing each 
of them by their respective harmonic order [i.e. the index number in Fig. [7|^c)]. In the case of 
the modes N, O, and P, the curves are perfectly on top of each other. For mode M, the scaling 
is slightly less good. Harmonics can be generated by several mechanisms. In the following, we 
will briefly discuss the parametric effect, electrical nonlinearities, and mechanical nonlinearities as 
possible origin of the harmonics. 

A common way to explain harmonics is based on the parametric effect [SI [7]. A mode at a 
resonance frequency /o can be actuated by varying the resonator spring constant /c at a frequency 
2/0/ j, where j is an integer > 1. It is easy to parametrically drive a nanotube resonator with a 
gate voltage, because /o (and therefore k) is widely tunable with Vg [8J. There exists a threshold 
above which the motion sets in: this threshold takes the form = {fodVg/dfo)/Q, where 
dfo/dVg is the change of the resonance frequency with gate voltage. For mode we get dfo/dVg = 
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FIG. 7: Map of resonance frequencies as a function of gate voltage measured with the FM 
technique, (a) Frequency modulation measurement at large driving force (obtained by measuring 
function of / and Vg with V^^ = 20 mV). Colour scale: (black) to 1 nA (dark red), (b) 
Low frequency range at even higher driving force {V^^ — 40 mV). Colour scale: (black) to 0.2 nA 
(dark red), (c) Map of all detected modes and harmonics. The number in the label designates the 
harmonic order of a resonance. Black lines correspond to resonances that cannot be assigned to a 
detected mode. 

36MHz/V close to = 4V. Together with the quality factor Q - 350 and /o - 124 MHz, 
this yields a threshold of 10 mV. This is consistent with the harmonics for mode N in Fig. [7|^a) 
^yac _ 20 mV). However, parametric excitation cannot account for the harmonics of order 4 and 
5 of mode M. 

A second scenario for harmonics is related to electrical nonlinearities in the circuit. Nonlinear- 
ities in current- volt age characteristics can generate forces at 2, 3, 4,... times the frequency of the 
applied V^^ and thus lead to harmonics with an index number n < 1. However, harmonics with 
an index number n > 1 are unlikely to have an electrical origin. 

Mechanical nonlinearities are predicted to give rise to harmonics. For instance, the quadratic 
nonlinear force can cause harmonics at 2/o and /o/2, and the cubic nonlinear force can induce 
harmonics at 3/o and /o/3 [6^ 9J. The combination of the quadratic and the cubic nonlinear forces 
can lead to harmonics with index 1/2, 1/3, 1/4, 2, 3, 2/3, 3/2,... [6J, which is in agreement with 
our measurements. 
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FIG. 8: Comparison of the modes with their harmonics. Each harmonic is divided by its respective 
order [corresponding to the index in Fig. [7|^c)]. (a) Mode M and harmonics. The subharmonics 
(M2, M4, and M5) deviate shghtly at low values of Vg. (b) Mode N and harmonics, (c) Mode O 
and harmonics, (d) Mode P and harmonics. 

In conclusion, the origin of the harmonics is not clear at the moment and this calls for future 
work. 



VII. SIMULATIONS 

We perform finite element simulations with ANSYS (R) Release 13.0 to reproduce the Vg de- 
pendence of the resonance frequencies. The mechanical properties of carbon nanotubes are well 
described by continuum elasticity and are independent of the chirality [lOj. For these simulations, 
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FIG. 9: Static shape of the nanotube. (a) Scanning electron micrograph of the nanotube. The 
electrodes are not voltage biased. Yellow dots mark the data points that we used to model the 
nanotube shape (for Vg = OV). Scale bar: 600 nm. (b) Points extracted from the electron micro- 
graph (black dots) and polynomial fit used for the nanotube model (red line), x is the coordinate 
along the axis connecting the two clamping points, y is the coordinate standing orthogonal to it. 

we use a tube with length L = 1.77 /xm, radius r = 1.5 nm, wall thickness Ar = 0.335 nm, mass 
density p = 2300 kg/m^, and Young modulus E = 1 TPa. We use the shape of the nanotube 
extracted from the scanning electron micrograph of the device (Fig. [9]). We assume that the static 
deformation is only in the horizontal plane when the device is not voltage biased. 

We use the 1-D BEAM188 element suitable for analyzing slender beam structures. 1-D 
BEAM 188 is a two-node element in 3-D and has six degrees of freedom at each node: transla- 
tions in the x, and z directions and rotations about the x, y, and z directions. A circular tube 
section is associated to the element by providing the inner and the outer radii. The constraints at 
the clamping points are fixed by setting all degrees of freedom to zero. Because of the high aspect 
ratio of nanotubes, we use point-like clamping conditions. The effect of the angle of the nanotube 
with respect to the electrodes is thus not accounted for. The electrostatic force induced by the 
dc voltage applied on the gate electrode, F^, is analytically calculated. For each gate voltage, 
the static deformation of the nanotube is calculated by performing a nonlinear structural analysis. 
This static solution is used as a base for the modal analysis. 
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VIII. EULER-BERNOULLI EQUATION: MODE FREQUENCIES 



In the following, we demonstrate that a good qualitative understanding of the motion of a sus- 
pended nanotube is possible from the Euler-Bernoulli equation. This approach provides analytical 
solutions that capture the behaviour of the system and that are similar to the solutions of the 
finite-element simulations performed with ANSYS. 

The Euler-Bernoulli equation for the static and dynamic displacement of a thin beam reads 

where p is the mass density, S the beam's cross-sectional area, z the displacement, t the time, E 
the Young modulus, / the second moment of inertia about the longitudinal axis, x the coordinate 
along the axis. To the built-in tension, L the resonator length, and g{t) a unit length force that 
accounts for the effect of the gate electrode in our experiment. We divide the displacement into a 
static and a dynamic component. 



dx'^ 



(4) 



z{x,t) = Zs(j)s{x) + Zi{t)(f)i{x) (5) 

where Zg is the maximum static displacement, zi is the maximum dynamic displacement, and 
05 (x), 01 (x) are the normalized static and dynamic profiles along the beam. 

In a first example, we develop Eq. [4] for the case of a single mechanical mode with 

05(x) = 01 (x) = sin(7rx/L) (6) 

where both the static and dynamic profiles are in the plane perpendicular to the gate electrode. 
This mode profile is strictly correct for negligible bending rigidity {EI 0). We choose this ansatz 
because it allows a simple analytical treatment of our problem. Moreover, we will see at the end of 
this section that this ansatz predicts gate voltage dependencies of the resonance frequencies that 
are in qualitative agreement with finite element simulations. We insert Eq. [5]and Eq. [6]into Eq. |4j 
multiply Eq. [4] by 0i(x), and integrate it from to L to get 
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In a static equilibrium position, the sum of the static terms in the first bracket on the right 
hand side of Eq. [7| is zero: 



^/..y +To..(-) --m. 



(8) 



The other terms of Eq. [7| can be rewritten in the usual form of a Newton equation of motion, 



(fzi{t) 



= -(jjlzi{t) - a2zl{t) - a^zl{t). 



(9) 



From a comparison of Eq. [7|to Eq. |9j we see that both 0^2 and 0^3 are positive, and that a2 (x Zg 
(meaning that a2 will vanish if the tube is straight). Further, both nonlinear coefficients are 
inversely proportional to and will become large for a short tube. They arise from additional 
tension that is generated when the beam bends (the integral term (^)^^^ [4] becomes 



nonzero . 



We repeat this calculation for other modes. The profile of the second mode, 02 (^), has the same 
shape as 0i(x), but stands orthogonal to it, moving parallel to the gate electrode. In this case, the 
final equation describing the second mode is somewhat simpler than that for the first mode: 



df2 
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zlit). (10) 



Again, 0^3 is positive, but here a2 = because the beam features no static bending in the 
direction of its vibrations. The linear restoring force of zi{t) is always larger than that of 2:2 (t), 
causing cji > 002- 

We assume that the third and fourth modes have the profile 



(l)s{x) = 04(x) = sin(27rx/L), 



(11) 



moving towards and parallel to the gate electrode, respectively. The resonance frequencies of 
the two modes are degenerate. Neither of these modes are detected in our experiment due to the 
antisymmetrical mode profile. The solution in this case reads (with i = 3 or 4) 
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We assume that the fifth and sixth modes are analogous to the third and fourth, but with 



05(x) = 06(^) — sin(37rx/L). 
Again, they are degenerate. We get (with j = 5 or 6) 



(13) 
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At this point, we can calculate the resonance frequencies of the first six modes as a function of 
Vg. For this, we determine the static displacement that provides an equilibrium of forces by solving 
Eq. Is) where 



(15) 



is the unit length force due to Vg^ and 



27r£o 



(16) 



d\n{2d/rY 

is the differentiation of the unit length capacitance with respect to the displacement. Here, d is 
the distance between the nanotube and the gate electrode, = 8-85 • 10~^^ Fm~^ is the electrical 
permittivity of free space, and r is the nanotube radius. We measure d — 370 nm by atomic force 
microscopy (AFM), but the large surface roughness of the electrodes of this device does not allow 
the measurement of r. We therefore chose a typical value from earlier nanotubes grown by the 
same method (r — 1.5nm). We also use E — 1 TPa, p — 2300 kg/m [Tl|, L — 1.77/im (measured 
by AFM), and wall thickness Ar = 0.335 nm. From these values, we calculate the second moment 
of inertia / = 2.928 • 10~^^kgm^ and the tubular cross-section S — 7T{{r + Ar)^ — (r — Ar)^) = 
3.157 • 10~^^ m^. The only free parameter, Tq, is chosen by optimizing the agreement of the mode 
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FIG. 10: Calculation of resonance frequencies as a function of Vg for a straight nanotube. (a) 
Resonance frequencies of the first six modes calculated from the Euler-Bernoulli equation, red: zi, 
green: black dashed: zs and Z4 (not detected in our experiment due to the antisymmetrical mode 
profile), blue: 2:5, orange: zq. Inset: shape of the nanotube when = V. (b) ANSYS simulation 
for the same set of parameters. The results are very similar to the analytical calculations. 

frequencies to the experimental results. We obtain To = 0.1 nN. The resulting Zg as a function of 
Vg is shown in Fig. 2(e) of the main text. Finite element calculations with a straight tube yield 
values of Zg close to the results of Eq.jsj see dashed line in Fig. 2(e) of the main text. The resonance 
frequencies of the six first modes are depicted in Fig. [To|(a). 



IX. EULER-BERNOULLI EQUATION: COUPLING BETWEEN MODES 

In the previous section, we have disregarded terms that couple different modes. The coupling 
has the same origin as the quadratic and cubic nonlinearities (0:2 and 0^3). It arises from the tension 
that is induced in a mode when another mode oscillates (through the integral term (^)^ in 
Eq. [4]). As a consequence, the coupling coefficients are of the same order of magnitude as a2 and 
0^3. Assuming that z{x^t) — Zs(l)s{x) + zi(j)i{x) + Z2(l)2{x)^ Eq. [4]leads to 

(fzA it) 

= -uj^ziit) - a2zl{t) - aszfit) - ^22Z^{t) - 6122^1 (t)z|(t) (17) 
where cjq, 0^2, and 0^3 are given by the expressions in Eq. jrj /322 = ^ (r)^ ^122 — ^ (r)^ 
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(all other summands /3jk and ej^i that appear in Eq. 1 of the main text are zero) . The equation of 
motion for the lowest mode moving parallel to the gate electrode is 



-UJoZ2it) - a3Z2it) - Pl2Zlit)z2it) - eii2Z2it) zf{t) 



(18) 



where uq and 0^3 are given by Eq. 10 and the couphng coefficients are /3i2 — j- (x) and 
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X. ONSET OF NONLINEARITY 

Nonhnear effects set in when zq reaches a critical value Zc- From Ref. [7J, we infer Zc = 
1.24:UJo/ ^/Q~\a\ assuming nonlinear damping to be negligible. Here, a is the effective nonlinear 
coefficient (see next section) that can be extracted from the backbone function connecting all 
resonance peaks at different driving amplitudes. From Ref. [7J, we have 

3azn 

where (jUmaxf^^ is the frequency where the amplitude is largest and (jJq/27v is the resonance 
frequency in the linear regime. In Fig. 1(d) of the main text, we have uJrnax — 27r • 46.85 MHz, 
(jjQ = 27r-46.35MHz, zq = 8.7nm, and thus get a = 3.2-10^-'^ m~^s~^. Using Q = 230 extracted from 
the resonance width in Fig. 1(c) of the main text, this leads to a critical amplitude of Zc = 4.2 nm, 
which is consistent with zq = 3.2 nm in Fig. 1(c) (where the resonance displays no hysteresis) and 
with Zq = 8.7 nm in Fig. 1(d) (where there is hysteresis). We repeat the same procedure for the 
data in Fig. 1(f) of the main text. Here, ujmax = 27r • 123.5 MHz, ujq = 27r • 124.25 MHz, zq = 2.1 nm, 
and we get a — —2.2 • 10^^m~^s~^. With Q — 354, we calculate Zc — 1.1 nm, which again is 
consistent with the results in Fig. 1(e) and (f) ( where zq = 0.9 nm and 2.1 nm, respectively). 



XI. STATIC DISPLACEMENT 

The reversal of the asymmetry of the resonance between Fig. 1(d) and (f) of the main text is 
due to a sign change of the effective nonlinearity a [3 [12] which depends on the quadratic and 
cubic coefficients in Eq. |9] as ^ 



10 _2 2 



(20) 



A comparison of Eq. [7] and Eq. [9] reveals that 



^2 = -;-z. 
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(21) 



and 



^3 



(22) 



E /7r\4 

The asymmetry of the resonance can be used to estimate the static displacement Zg of the 
resonator at different values of Vg. With a known from the estimations in the last section, we can 



insert Eq. 21 and Eq. 22 into Eq. 20 in order to obtain Zg. With E = 1 TPa, p = 2300 kgm ^, and 
L — 1.77 /xm, we get Zg — 2.8 nm for Vg — 1.5 V, and Zg — 13 nm for = 4V. 



XII. ELECTROSTATIC NONLINEARITIES 

In the previous sections, we have shown that the sign change of the effective nonlinearity a 
(from a positive to a negative value) as a function of Vg is consistent with the expected increase 
of the quadratic nonlinearity 0^2 • A negative a could also have an electrostatic origin [7^ T2]. 
The quadratic and cubic nonlinearities are af = ~^^'''^g ^3 ~ ~^^""^g ^ where m 

is the effective mass of the resonator and C" and C"" are the third and fourth derivatives of 
the capacitance with respect to displacement. The values calculated for = 4V are = 
—4.4 • 10^^ m~^s~^ and — —4.0 • 10^^ m~^s~^, which are at least 3 orders of magnitude smaller 
than the mechanical nonlinearities above. We therefore neglect electrostatic nonlinearities in the 
analysis of our experiment. 



XIII. RESONANCE LINESHAPES 



In Fig. 4 of the main manuscript, we show resonance lineshapes with exotic features that we 
associate with mechanical coupling between commensurate modes of the nanotube. Here, we show 
how these features vanish and how a usual Duffing nonlinearity is recovered when we reduce the 




FIG. 11: Resonance lineshape for different driving forces. The resonance is measured with the 
two-source method for Vg = 1.88V and = 0.28mV. = 0.56mV in (a), = 2.2mV in 
(b), V^^ — 3.4 mV in (c), and V^^ — 4.5 mV in (d). A conventional Duffing nonhnearity (with 
q; < 0) is recovered at the lowest driving force. Sweeps are performed with decreasing frequency. 
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FIG. 12: Mixing current as a function of / and Vg for different driving forces. I mix 

is measured 

with the two-source technique with V"/^ = 0.28 mV. V"^^ = 0.56 mV in (a), V"^^ = 1.1 mV in (b), 
V^^ = 2.2 mV in (c), V^^ = 3.4mV in (d), and V^^ = 4.5 mV in (e). Sweeps are performed with 
decreasing frequency. Panel (b) corresponds to Fig. 4(f) of the main text. 
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FIG. 13: Increase of the resonance width with V^^. We measure the resonance of mode N at 
Vg = 1.5V with the FM technique for = O.lmV (black), 1/^^ = 0.12mV (red), 1/^^ = O.lSmV 
(blue), and V^^ = 0.2 mV (green). The distance between the two minima flanking the resonance 
peak (solid bars) corresponds to the resonance width A/ = fo/Q- 



driving force. The different panels in Fig. [TT] correspond to the same resonance for different driving 
forces (oc V^^). 

A second example with increasing driving force is displayed in Fig. [T2j Here, we can observe the 
evolution from an almost unbroken line (resonance frequency versus Vg) in Fig. [T2|(a) to a highly 
exotic response with multiple peaks and dips as a function of / in Fig. fl2Fe). 



XIV. NONLINEAR DAMPING 

We recently reported that the resonance width A/ = fo/Q of nanotube and graphene mechan- 
ical resonators can depend on the driving voltage V^^. We attributed this phenomenon to the 
nonlinear damping force r]z^z [13j. In the present device, the dominant bistability behaviour pre- 
vents observing nonlinear damping above V^^ ^ 0.2 mV. Nontheless, we found an increase of the 



resonance width for V^^ < 0.2 mV (Fig. 13). In this measurement, we use the frequency modulation 



technique, which produces two characteristic minima flanking the resonance peak. The separation 
of these minima corresponds to A/. A clear increase of A/ is seen between V^^ = O.lmV (below 
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FIG. 14: Map of resonance frequencies as a function of V^. (a) Frequency modulation measurement 
at low driving force (obtained by measuring function of / and Vg with V^^ = 4mV). Three 

resonances are detected with a large signal, a fourth one shows up faintly. Colour scale: (black) to 
0.1 nA (dark red), (b) Same measurement with a larger frequency range and V^^ = 40 mV. Colour 
scale: (black) to 0.1 nA (dark red), (c) Schematic of the map of the resonance frequencies as a 
function of Vg. Resonances that cannot be assigned to a mode are drawn in black. Black arrows 
mark regions where modes O and P are commensurate and the resonances lineshapes become 
exotic. 




4.2 4.4 4.6 

Vg (V) 

FIG. 15: Modal coupling in the second device. Measurement with the two-source method (ob- 
tained by measuring Imix ^is a function of / and Vg with V^^ — 2.8 mV and V^^ — 0.28 mV). Colour 
scale: (black) to 0.1 nA (dark red). 
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this driving voltage the signal vanishes) and V^^ = 0.2 mV. 



XV. ADDITIONAL DEVICE 



A second nanotube resonator exhibits similar behaviour as the one discussed so far. The results 



of the second device are summarized in Fig. 14 and Fig. [TSj Three modes are clearly visible at low 
driving force, while a fourth resonance shows up faintly [Fig. [T4|(a)]. In these measurements, we 
use the FM technique with V^^ = 4mV. With a larger driving force (V^^ = 40 mV), many more 
resonances appear [Fig. [l4]^b)]. With this device, we have mapped the frequency spectrum up to 
500 MHz, and all detected resonances are depicted in Fig. [T4|(c) and labelled according to the most 
probable harmonic spectrum. We can identify regions where two modes are commensurate or nearly 
commensurate and the resonance lineshapes become exotic (black arrows). Around Vg = 4.4V, 
mode P has exactly twice the frequency of mode O (Fig. [l5]). There, we observe a discontinuity 
in the map of the resonance frequency as a function of Vg. 



[1] J. Cao, Q. Wang, and H. Dai, Nature Mat. 4, 745 (2005). 

[2] A. K. Hiittel, G. A. Steele, B. Witkamp, M. Foot, L. P. Kouwenhoven, and H. S. J. van der 

Zant, Nano Lett. 9, 2547 (2009). 
[3] G. A. Steele, A. K. Hiittel, B. Witkamp, M. Foot, H. B. Meerwaldt, L. P. Kouwenhoven, and 

H. S. J. van der Zant, Science 325, 1103 (2009). 
[4] V. Gouttenoire, T. Barois, S. Perisanu, J.-L. Leclercq, S. T. Purcell, P. Vincent, and A. Ayari, 

Small 6, 1060 (2010). 

[5] V. Sazonova, Y. Yaish, H. Ustiinel, D. Roundy, T. A. Arias, and P. L. McEuen, Nature 

(London) 431, 284 (2004). 
[6] A. H. Nayfeh and D. T. Mook, Nonlinear Oscillations (Wiley- VCH, New York, 1979). 
[7] R. Lifshitz and M. C. Cross, Reviews of Nonlinear Dynamics and Complexity 1 (Wiley- VCH, 

New York, 2008), available at www.tau.ac.il/ ronlif/pubs/RNDCl-l-2008-preprint.pdf. 
[8] A. Eichler, J. Chaste, J., Moser, and A. Bachtold, Nano Lett. 11, 2699 (2011). 
[9] M. I. Dykman, D. G. Luchinsky, R. Mannella, P. V. E. McClintock, S. M. Soskin, and N. D. 

Stein, Phys. Rev. E 54, 2366 (1996). 



[10] K. N. Kudin, G. E. Scuseria, and B. 1. Yakobson, Phys. Rev. B 64, 235406 (2001). 
[11] C. Lee, X. Wei, J. W. Kysar, and J. Hone, Science 321, 385 (2008). 

[12] I. Kozinsky, H. W. Ch. Postma, I. Bargatin, and M. L. Roukes, Appl. Phys. Lett. 88, 253101 
(2006). 

[13] A. Eichler, J. Moser, J. Chaste, M. Zdrojek, L Wilson-Rae, and A. Bachtold, Nature Nan- 
otechnol. 6, 339 (2011). 



